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ABSTRACT 

Stars and dark matter account for most of the mass of early- type galaxies, but uncertainties in 
the stellar population and the dark matter profile make it challenging to distinguish between the two 
components. Nevertheless, precise observations of stellar and dark matter are extremely valuable for 
testing the many models of structure formation and evolution. We present a measurement of the stellar 
mass and inner slope of the dark matter halo of a massive early- type galaxy at z = 0.222. The galaxy 
is the foreground deflector of the double Einstein ring gravitational lens system SDSSJ0946+1006, 
also known as the "Jackpot" . By combining the tools of lensing and dynamics we first constrain the 
mean slope of the total mass density profile (ptot oc ) within the radius of the outer ring to be 
7' = 1.98 ±0.02 ±0.01. Then we obtain a bulge-halo decomposition, assuming a power-law form for the 
dark matter halo. Our analysis yields 7dm = 1.7 ± 0.2 for the inner slope of the dark matter profile, 
in agreement with theoretical findings on the distribution of dark matter in ellipticals, and a stellar 
mass from lensing and dynamics = 5.5^5 3 ^ lO^^M©. By comparing this measurement with 

stellar masses inferred from stellar population synthesis fitting we find that a Salpeter IMF provides 
a good description of the stellar population of the lens while the probability of the IMF being heavier 
than Chabrier is 95%. Our data suggest that growth by accretion of small systems from a compact 
red nugget is a plausible formation scenario for this object. 

Subject headings: Galaxies: early-type — Galaxies: structure — dark matter — gravitational lensing 



1. INTRODUCTION 

What are the main physical processes that shape early- 
type galaxies (ETGs)? Understanding the formation and 
evolution of early-type galaxies is a fundamental piece in 
the cosmo logical puzzle. Any model that aims at pro- 
viding a description of the Universe as a whole must be 
able to reproduce the observed characteristics of these 
objects. 

ETGs are observed to hold tig ht scaling relations , 
such as the Fundamental Plane (Dressier et al.l 119871 : 
iDjorgovski fc Davis 1987) and the correlation between 
the mass of the cen tral black hole and global galac- 
tic properties (.Fer rarese fc Merrittll2000l : iGebhardt et al.l 
120001 : iMarconifc Hunt 2^. The inner density pro- 
file of their total mass is measured to be very close 
to isothermal, in th e so-call ed ^^bulge-halo conspiracy" 
(iTreu fc Koopmans' ^04; Ko opmans et al.l [2QQ6L 12001 
lAuger et al. 2010a). Moreover, they appear to undergo 
a significant evolution in size, from being very compact 
in the early {z ^ 2) Universe to the more diffuse objects 
that we observe at more recent times (|van der Wei et all 
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l2QQ8l : lvan Dokkum et ani2QQ8l : iNewman et al.llMol ). Fi- 
nally, studies of massive ETGs seem to favor a heav- 
ier stellar initial mass function (IMF: iGrillo et al.1 



Treu et all 120101: Ivan Dokkum fc Conrovl [2QTqI 



2009; 



2011 



Auger et al .' '2010bl: ISDiniello et al l 120 111) than for spi- 
ral galaxies (Beh fc de Jond 120011 : iDutton et all 120111 : 
ISuvu et al.ll201 1). 

These characteristics should reflect a common mecha- 
nism that drives ETGs towards the tight relations ob- 
served at low redshifts {z < 0.3) during their forma- 
tion and evolution. Mergers with other galaxies are 
likely to be one of the key processes in the history 
of ETGs. Mergers are believe d to be at the basis 
of th e formation (Kormendv fc Sanderslll992l : IHernquistI 
119921 : fShier fc Fischer 1998). to be involved in the black 



hole scaling relation ( Haehnelt fc Kauffmannll2000l: iPeng 
2007f) and to drive the size evolution (jHopkins et al ' 
2009al : lOser et al.l I2012D of ETGs. However, the pic- 



ture is compHcated by the many physical processes that 
are present during the evolution of galaxies, such as gas 
cooling, feedback from stars and AGN. Moreover, an 
important fraction of the mas s of ETGs i s accounted 
for b y dark matter (Bertin et al.l IT991: iFranx et al.l 
|l99i; "Gerhard et al." 2001,; >Treu fc Ko opmansl [2QQ1 
iBarnabe et al. 2011) whose nature is still unknown. Un- 
derstanding the interplay of baryonic and dark matter 
and how they act to produce the observed structural 
characteristics is essential to comprehend the evolution 
of ETGs, but is today a challenging task. 

Theoretical models and simulations with a variety of 
physical ingredients have been set up to try to repro- 
duce the obse r yables of ETGs (iGustafsson et all 120061: 
iHopkins et all 120101 : iSchave et al.l 120101 : iDuffv et al.l 



2 



Sonnenfeld et al. 



l2QlQf ). Although simulations seem to be able to capture 
the general characteristic s of ETGs (e.g. iHopkins et al.l 
l2QQ9bl : iCiotti et al]l2QlQ[ ). quantitatively matching the 
entire set of observables proved to be difficult, often 
requiring an ad h o c tuning of the mo del parameters 
jNipoti et~aI ||2QQ9l: I Hopkins et aIll2QlQ[ ). For example, 
iDuffv et ldl 1^10) explored a variety of aspects of bary- 
onic physics such as gas cooling, feedback from stars and 
AGN, finding that on the one hand the observed inner 
slopes of massive ETGs are reproduced if the feedback is 
weak, but on the other hand a strong feedback is needed 
to match the measured stellar masses. 

Improving the quality of the observation of ETGs and 
introducing more constraints can help us to discriminate 
between the wealth of currently viable scenarios for their 
history. Two characteristics of ETGs in particular are 
still not known with sufficient precision and leave room 
for significant improvement in their observational deter- 
mination: the stellar mass and the density profile of the 
dark matter halo. The stellar mass is degenerate with 
the stellar IMF with respect to constraints from the in- 
tegrated light distribution and colors. Breaking this de- 
generacy can help determining the star formation history 
and the content of the baryonic mass in ETGs. The den- 
sity profile of the dark matter distribution is sensitive to 
the physical processes that take place during the forma- 
tion and evolution of ETGs. Therefore, measuring the 
profile of dark matter halos in ETGs is a powerful means 
for testing the various theoretical models. 

Part of the difficulty in comparing simulations with 
observations are, of course, due to the fact that dark 
matter, which accounts for a significant fraction of the 
mass of a typical galaxy, is not directly observable. Grav- 
itational lensing is, in this aspect, a very powerful tool, 
being sensitive to the gravitational pull of matter inde- 
pendently on its interaction with light (e.g.. lTreull2QlQL 
and references therein). 

Lensing surveys indeed played a crucial role in uncov- 
ering physical characterics of early- type galaxies, such 
as theii^ average density profile and dark matt er frac- 
tion (Koopmans et al.l I2QQ6L I2QQ9I : |Auger et al| l2QlQal : 
iBarnabe et al. 20111) or t he IMF of their stellar popu- 
lation (jGrillo et al.l I2QQ9I : iTreu et al.|[2QlQl : lAuger et al.] 
l2QlQb[ ). Measurements of the inner slope of a dark mat- 
ter h alo have so far been obtained for a few c luster lenses 
(e.g. "Sand et al.l I2QQ8I : iNewman etHI l2QQ9f ). for which 
constraints from multiple lensed sources are available. 
For typical early- type galaxy strong lenses however, there 
are residual degeneracies between anistropy, stellar mass 
to light ratio and inner slope of the d ark matter halo and 
therefore the constraints are weak (i Koopmans fc Treul 
[2003; Treu & Koopmans 2004). For this reason, previous 
studies have adopted theoretically mot ivated mass den- 
sity profiles for the dark matter halo (|Treu et al.ll2010l : 
lAuger et al.ll201Qbl ). rather than free power laws. 

Here we present a detailed study of an early-type 
galaxy at redshift z — 0.222. The galaxy is the strong 
gravitational lens of the system SDSSJ0946+1006, part 
of the SLAGS sample (Bolton et al. 2004). This ETG 
is special in that it lenses two sources at different red- 
shifts, creating two nearly complete Einstein rings of 
different radii. For this reason, the system is also re- 
ferred to as the "Jackpot". The first lensed source is 
at redshift Zs\ = 0.609, while there is no spectroscopic 



measurement of the redshift of the second ring. Thanks 
to the presence of the two rings, this system provides 
more information than typical gravitational lenses, de- 
spite the lack of the second source redshift . A first study 
of SDSSJ0946+1006 was carried out by Gavazzi e t al.l 
(2008, Paper VI). An indepe ndent lensing analys is of this 
system was performed by iVegetti et al.l (|2010f ). which 
led to the discovery of a small satellite with no visible 
counterpart. Here we include new high-quality photom- 
etry obtained with the Hubble Space Telescope (hereafter 
HST) and new deep and spatially resolved spectroscopy 
obtained at the Keck Telescope. The goal of our study 
is to separate the contribution of dark and stellar matter 
to the total mass of the lens, making as few assumptions 
as possible about the density profile of the dark matter 
halo. This task is achieved by combining lensing and 
dynamics information. Unlike typical early-type galaxy 
lenses, the wealth of information provided by this system 
allows us to determine both the mass of the stellar bulge 
and the inner slope of the dark matter halo. Thanks to 
the multi-band HST photometry we are able to obtain a 
photometric redshift of the outer ring, that is necessary 
for improving the constraints from the lensing data, and 
to infer stellar masses from stellar population synthesis 
(SPS) fitting. The comparison between this measure- 
ment of the stellar mass and the one obtained through 
lensing and dynamics allows us to constrain the IMF of 
the stars in the lens. This is the most robust measure- 
ment of the inner slope of the dark matter halo and IMF 
of an isolated massive ETG. 

The structure of this paper is the following. In Sections 
2 and 3 we describe the new photometric and spectro- 
scopic data, respectively. Our measurement of the pho- 
tometric redshift of the outer ring is presented in Section 
4. In Section 5 we describe measurements of the stellar 
mass of the lens from stellar population synthesis fitting. 
Section 6 describes a lensing and dynamics model as- 
suming a power-law density profile for the total density 
profile of the lens. In Section 7 we present the bulge- 
halo decomposition of the lens. We discuss our results in 
Section 8 and summarize in Section 9. 

Throughout the paper we assume the following 
values for the cosmological parameters: Hq = 
70 km s"^Mpc"\ = 0.3, = 0.7. Magnitudes are 
expressed in the AB system, images are North-up and 
position angles are in degrees East of North. In show- 
ing our results we display posterior PDFs in multiple 
projections wherever possible, but when giving a point 
estimate of an inferred parameter we quote the position 
of the peak of its one-dimensional marginalised distri- 
bution, with uncertainties defined by the 68% credible 
region. 

2. MULTICOLOR HST PHOTOMETRY 

We present HST images of the lens system 
SDSSJ0946+1006 in five different bands. In Paper VI 
we reported results based on an ACS F814W image only. 
Images in WFPC2 F606W and WFC3-IR F160W (Cy- 
cle 16, Program 11202, PI K oopmans) were available for 
Paper IX (Aug er et al.ll2009l ). In addition to those data, 
we now have WFC3 images in F438W and F336W bands 
(Cycle 17, Program 11701, PI Treu). Table [T] summarizes 
the observations. This section describes the data reduc- 
tion process (§ 12. ip and the photometric properties we 
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Fig. 1. — Gravitational lens system SDSSJ0946+1006 in a com- 
bination of F814W, F606W and F336W HST images. 



TABLE 1 
Summary of the HST observations. 



Instrument 


Filter 


Exp. time 


A^ exp 


Date 


WFC3 IR 


F160W 


2397 s 


4 


09/12/2009 


ACS 


F814W 


2096 s 


4 


3/11/2006 


WFPC2 


F606W 


4400 s 


4 


18/12/2009 


WFC3 UVIS 


F438W 


2520 s 


4 


20/03/2010 


WFC3 UVIS 


F336W 


5772 s 


4 


20/03/2010 



derived for the lens galaxy (§ \2.2\ . For conciseness, we 
sometimes refer to the F160W, F814W, F606W, F438W, 
F336W bands as H, I, V, B, U respectively. A color com- 
posite image of the lens system is shown in Figure [T] 

2.1. Data reduction 

The data are treated with the standard HST reduction 
pipeline. For each image, frames are coadded and re- 
sampled in a uniform pixel scale using the software MUL- 
TIDRIZZLE CFruc hter et al.ll2QQ9l ). Pixel sizes are 0.10'^ 
for the F160W image, 0.050'' for the F814W and F606W 
images, and 0.0396'' for the F438W and F336W im- 
ages. The images are then brought to the same orienta- 
tion and 0.050" p ixel scale by using the software SWARP 
(|Bertin efall[2002 ). The PSF of each image is estimated 
from stars in the field. 

2.2. Lens galaxy properties 

The brightness distribution of the main lens galaxy is 
first obtained by fitting Sersic profiles to the data. This 
task is achieved with the software SPASMOID, develope d 
by M. W. Auger and described by Bennert et al. (2011). 
Spasmoid fits the data in all the bands simultaneously 
with a unique model, determining total magnitude and 
colors of the galaxy at once. By using a single Sersic com- 
ponent we find a best-fit profile described by a Sersic 
index n = 6.0, axis ratio q = 0.95 and effective radius 
reff = 2.93". However, the residuals left by this single- 



component fit are rather large. Consequently, we add a 
second component, allowing for the position angle of the 
major axes of the two profiles to be different but impos- 
ing a common centroid. In the fitting process, the light 
from the rings is masked out manually. This procedure 
gives robust estimates of the colors of the lens, rather in- 
dependent from the model adopted to describe the data. 
Color information will be used in Section 5 to constrain 
the stellar population. In Fig. [2] we show the images of 
the system in the five bands, before and after subtracting 
the main lens. Residuals are on the order of a few percent 
in the F814W band image. Table [2] reports the best-fit 
structural parameters of the model, while the best-fit col- 
ors are given in Table |4j It is worth pointing out that 
the major axes of the two components are almost per- 
pendicular, and that the mean surface brightness within 
the effective radius of component 1 is a factor ~ 30 larger 
than that of component 2. The measured magnitude in 
the F814W band is c onsistent with the value reported by 
iGavazzi et al.l (|2008f ) for the same object. 

In order to both explore model-dependent systematic 
errors and obtain a computationally more tractable de- 
scription of the light profile for our lensing analysis, 
we also model the lens light with the following surface 
brightness distribution: 



I{x,y) Ic 



1 



1 



(1) 



where R^ = x^/g + qy^. This profile corresponds to a 
truncated pseudoisothermal elliptical mass distribution 
(tPIEMD) in 3d, with Tc and rt corresponding to the 
core radius and truncation radius respectively. Note that 
the number of parameters of the model is the same as 
that of a Sersic profile. Two components are used, as in 
the Sersic case. The best-fit parameters are reported in 
Table O Both the double-Sersic and the double-tPIEMD 
profiles fit well the photometry of the lens, with residuals 
within the outer ring on the order of a few percent in the 
F814W band (see Figure ED- 

The inferred total magnitude in the two models is dif- 
ferent, but this is due to the different behavior at large 
radii, where there are no data. In fact, the magnitude 
within the inner ring is the same for the two models to 
within 0.01 mags and the inferred colors are consistent 
within the errors with those reported in Table HI 

The infrared F160W data reveal distorsions in the 
shape of the light distribution at large radii (see Figure 
[3|), a possible si gnature of tida l inter actions. As pre- 
viously noted bv IGavazzi et al.l ()2008f ). a galaxy in the 
neighborhood of the lens also shows signs of a tidal inter- 
action (see Figure [3]). It is possible that the two galaxies 
are undergoing a merger. This deviation from a regular 
light profile is located far from the probed by our lens- 
ing and dynamics measurements and is therefore not a 
concern for the accuracy of our models. The central part 
appears smooth to the few percent level and it is unlikely 
that the ongoing interaction would have an effect on its 
structure, given its deep potential well. However, as we 
will discuss in § 18.21 this feature provides an interesting 
clue to the formation mechanism of this galaxy. 

Another interesting feature is revealed by the image in 
the F336W (U) band, as there seems to be some struc- 
ture in the center of the lens (see Figure [2]). The fact that 
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Fig. 2. From top to bottom: HST F160W, F814W, F606W, F438W and F336W images of the lens system SDSSJ0946+1006 before 
(left column) and after (middle and right column) light subtraction. Middle column: light distribution modeled as a double Sersic profile, 
with parameters given in Table [2] Right column: light distribution modeled as a double tPIEMD profile, with parameters given in Table [3] 



this feature is clearly visible only in the U band, where 
the lens is fainter, may suggest that it is in fact a bluer 
object distinct from the central galaxy, or blue emission 
from an active nucleus. Alternatively, the observed de- 
tail could be the result of the presence of a dust lane 
that separates the light of the lens into two components 
at shorter wavelengths. In principle it could also be an 
additional image of the lensed sources. 

One way to discriminate between a blue object or a 
dust lane is to study the position of the centroid of the 
lens in the different bands. A blue object would shift 
the centroid towards itself at bluer wavelengths, while a 
dust lane would remove blue light, causing the appar- 



ent centroid to move away from it. When fitting for the 
centroid of the lens, this latter case is observed: the cen- 
troid moves by about one pixel towards the S-E in the 
F336W and F438W bands with respect to the F814W 
band. This is a significant effect given the subpixel ac- 
curacy of centroniding, and it suggests that dust is most 
likely the cause of the observed feature in the F336W 
band. A more detailed discussion of the dust issue is 
given in Appendix A. 



3. PHOTOMETRIC REDSHIFT OF THE OUTER RING 
3.1. Colors of the ring 
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TABLE 2 

Lens light distribution: double Sersic model 



Component 


^F814W 


^eff 


n 


q 


PA 


(^^)e,F814W 




(mag) 


(arcsec) 






(degrees) 


(mag arcsec"^) 


1 


18.38 ±0.20 


0.50 ±0.10 


2.34 ±0.50 


0.79 ±0.10 


63.0 ± 1.0 


18.87 ±0.10 


2 


17.44 ±0.10 


4.46 ±0.50 


1.60 ±0.50 


0.64 ±0.10 


-23.4 ± 1.0 


22.68 ±0.20 



Best fit parameters for the double-Sersic model surface brightness profile of the main lens: magnitude 
in the F814W band, effective radius, Sersic index (n), axis ratio (g), position angle of the major axis 
(East of North), effective surface brightness. Each line refers to one of the Sersic components of the 
model. The errors represent the typical range of values for the parameters allowed by the model. These 
errors are correlated: for example, an increase in the value of the Sersic index n results in a change of 
the effective radius to fit the observed slope in surface brightness. 



TABLE 3 

Lens light distribution: double tPIEMD model 



Component 


^F814W 


rc 


rt 




PA 


(SB)e,F814W 




(mag) 


(arcsec) 


(arcsec) 




(degrees) 


(mag arcsec ~-^) 


1 


18.75 ±0.20 


0.066 ±0.010 


0.50 ±0.05 


0.66 ±0.10 


63.0 ± 1.0 


19.24 ±0.10 


2 


17.15 ±0.10 


0.082 ±0.010 


6.05 ±0.10 


0.71 ±0.10 


-24.3 ±1.0 


22.48 ±0.20 



Best fit parameters for the double-tPIEMD model surface brightness profile of the main lens: magnitude 
in the F814W band, core radius (rc), truncation radius (rt), axis ratio (g), position angle of the major 
axis (East of North), effective surface brightness. 



Oolor 
I- H 
V - I 
B - V 
U - B 



TABLE 4 
Colors of the lens galaxy. 



Component 1 
1.16 ±0.05 
0.81 ±0.05 
2.36 ±0.20 
2.30 ±0.30 



Component 2 



0.86 ±0.05 
0.96 ±0.05 
1.52 ±0.05 
1.32 ±0.10 



dlobal 
0.96 ±0.05 
0.91 ±0.05 
1.73 ±0.05 
1.44 ±0.10 



TABLE 5 

AB COLORS OF THE OUTER RING 

1 - H 0.61 ±0.10 

V-I 0.21 ±0.10 

B-V 0.15 ±0.10 

U-B 0.53 ±0.10 




Fig. 3. — WFC3-IR F160W image of the lens system and its sur- 
roundings. Note the irregular shape of the faint stellar component 
at the outskirts of the lens galaxy (top of the image). At the bot- 
tom, a neighbor also shows signs of tidal disruption. Both these 
features may be the result of a close encounter between the two 
objects. 



One of the main goals of this study is to constrain 
better the mass distribution in the lens galaxy by ob- 
taining a photo-z of the outer ring. This task requires 
a measurement of the colors of the ring. A color map 
of the outer ring is obtained as follows. For each pair 
of neighboring bands, Ai, A2, we align the correspond- 
ing images and then convolve each image with the PSF 
of the neighboring band. In this way we obtain pairs of 
images with the same effective PSF, necessary to get an 
unbiased estimate of the color for each pixel. Global col- 
ors are then measured in the following way. For a given 
pair of bands, we select individual pixels with flux larger 
than the background by more than two sigma in both of 
the bands considered. We make the assumption that the 
source has spatially uniform colors and estimate them 
statistically by taking a weighted mean of the individual 
pixel colors. The measured values of the colors, corrected 
for galactic extinction, are reported in Table [H 

3.2. Measuring the photo-z 

To estimate the photometric redshift of the outer ring 
we make use of the software BPZ (Bayesian Photo-z; 
lBenftezll2QQQf ). Photo-z analysis consists of fitting syn- 
thetic SEDs to the observed colors. BPZ works in a 
Bayesian framework that allows us to combine the infer- 
ence with that from other pieces of information: given 
a prior probability distribution for the source redshift 
and galaxy type, BPZ calculates the probability of the 
source being at redshift Zs2 given its colors C and mag- 
nitude m, P(2;s2|C,m). The stellar templates used for 
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Fig. 4. — Solid line: posterior probability distribution function 
of the source redshift, as calculated with BPZ, assuming a prior on 
Zs2 from Hubble Deep Field North number counts. Overplotted 
are the levels corresponding to 68% and 95% enclosed probability. 
Dotted line: posterior PDF assuming a flat prior on Zs2- 

the SED fitting are described by 'Coe et al J ()20Q6r). Th e 
F814W magnitude is taken from Gavazzi et alJ (|2008l ), 
where the brightness distribution of the source was re- 
constructed after a lens modeling. The value adopted is 
therefore mF8i4w = 27.01 ± 0.19 

For the redshift distribution we u se a prior 
P(z|mF8i4w) suggested by Benitei (|2QQQf ) and based 
on number counts from the Hubble Deep Field North 
(HDFN). Figure m shows the redshift posterior probabil- 
ity distribution function P(2;52|C, mF8i4w)- The most 
likely redshift with 68% confidence interval is Zs2 = 
2.411q"2i. As will be shown later, this information is suf- 
ficient to put interesting constraints on the model of the 
lens system. We also calculated the photo-z assuming a 
fiat prior on Zs2^ and found a nearly identical result. Col- 
ors calculated with a different lens light subtraction, the 
double-Sersic model, yield the same photo-z well within 
the quoted uncertainties. 

4. KECK SPECTROSCOPY 

The data were collected during the nights of 2006 De- 
cember 23 and 24 with the LRIS instrument at the Keck 
Telescope 1. The original goals of the observations were 
to measure a velocity dispersion profile of the foreground 
defiector and to measure the redshift of the outer ring. 
The first goal was succesfully achieved, while we were 
not able to detect any spectroscopic signature from the 
farthest source. 

Because of the dual scope of our study, two different in- 
strumental setups were used. The first setup, used during 
the first night, was optimized for a better measurement 
of the velocity dispersion of the defiector. The wave- 
length range in the red detector, the one used for the 
measurement of a, was ~ 5700 — 7600A, bracketing im- 
portant absorption features in the rest frame of the lens 
dit z = 0.222. During the second night we centered the 
slit on the longest arc of the outer ring, and used a setup 
with a broader wavelength range, up to ~ 86OOA. A 
summary of the observations, with specifications on the 
setups used, is provided in Table [6l The 
spectrum of the system is shown in Fig. [5l There is no 



evidence for the presence of emission lines from objects 
other than the foreground lens and the inner ring. Given 
our measurement of the photo-z of the outer ring, we 
would expect Ly-a emission to fall around ~ 4150A, but 
it cannot be identified in our spectrum. We can put an 
upper limit of ~ 5 x 10~^^ erg cm~^ s~^ to the fiux in 
Ly-a from the source. 

4.1. Velocity dispersion 

The velocity dispersion of the main lens is measured 
by fitting stellar templates convolved with a Gaussian 
velocity distribution to the observed spectrum. This op- 
eration is carried out with a Monte Carlo Markov Chain 
approach, usi ng a code d e velop ed by M. W. Auger, and 
described by iSuvu et'ol] (j2010l ) . The rest frame wave- 
length range used for the fit is 5100 — 5850 A. For the stel- 
lar templates we used linear combinations of nine spectra 
from the INDO-US library, corresponding to K,G,F and 
A stars. The most prominent absorption feature in the 
wavelength range considered is Mgb (5175A). However, 
we experienced difficulty in finding a good fit to both 
Mgb and the rest of the spectrum. It is known that 
some galaxies have enhanced magnesium features in the 
spectrum that are not well repr oduced in standard stel- 
lar templates (jBarth et al.ll200^ . For this reason we de- 
cided to mask the Mgb absorption line out of the fitted 
spectrum. 

With the aim of obtaining a velocity dispersion profile, 
we measured <j in a set of apertures. The spatial position 
of the apertures was determined by fitting the centroid 
of the trace of the lens in the twodimensional spectra 
and it is accurate to ~ 0.02''. In Tableland Fig. [6] we 
report the measured values of a and of the mean velocity 
in each aperture, while in Fig. Owe show the fit in the 
central 0.42'' as an example. There is evidence for some 
rotation, with v^^^ <C cr^. 

5. STELLAR MASSES 

Here we present a measurement of the stellar mass of 
the foreground lens galaxy. The procedure adopted is 
the following: we fit stellar population synthesis mod- 
els to the observed spectral energy distribution (SED) 
of the galaxy. A measurement of this kind was al ready 
performed bv lGrillo et all §009) and I Auger erall (j2009L 
Paper IX) for the s ame ob ject. Their results agree within 
the errors. Grillo et all (|2009) used SDSS multiband 
photometry (li, ^, r, z, z bands) as their observed SED. In 
Paper IX, high resolution HST d ata was used, but on ly in 
two bands (F814W and F606W). lAuger et all (j2009[ ) also 
introduced a powerful statistical analysis method, based 
on Bayesian statistics that allows for physically mean- 
ingful priors on the model parameters as well as a full 
exploration of uncertainties and correlation between the 
inferred parameters. With five band HST photometry 
we can now extend the analysis of Paper IX, to obtain a 
more robust estimate of the stellar mass. 

The fitting method is the same as that developed by 
lAuger et al.l (|2009V and can be summarized as follows. 
Composite stellar p opu lation models are created from 
iBruzual fc Charlotl (|2Q0 3) stellar templates. The star 
formation history is modeled with a single exponentially 
decaying burst. The parameters of the model are age, 
metallicity, exponential burst timescale, dust reddening 
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TABLE 6 

Spectroscopic observations: summary 



Date 


Exp. time 


Slit width Dichroic 


Blue grism 


Red grating 


Red Ac 


Weather 


Seeing 


12/23/2006 


16200 


1.0" 560 


600/4000 


831/8200 


6819A 


Good 


0.8^^ 


12/24/2006 


12600 


1.0" 680 


300/5000 


831/8200 


7886A 


Good 


0.8^^ 



H/3 

z=0.609 




4000 



5000 6000 ^ 7000 

Observed wavelength (A) 



8000 



Fig. 5. — LRIS spectra of the Jackpot. Blue: data from the first night, 
from rectangular apertures 1" x 3.36^^ Dotted line: noise level. 

TABLE 7 
Velocity profile measurements 



Red: data from the second night. The two spectra are extracted 



Slit offset 




cr 


(arcsec) 


(km s-l) 


(km s-1) 


-1.05 


101 ± 21 


252 ± 25 


-0.84 


85 ±16 


273 ± 18 


-0.63 


62 ±11 


263 ± 14 


-0.42 


30 ±10 


278 ± 12 


-0.21 


20 ±10 


287 ± 11 


0.00 


0±9 


287 ± 11 


0.21 


-22± 11 


286 ± 11 


0.42 


-55± 12 


299 ± 13 


0.63 


-67± 13 


274 ± 15 


0.84 


-63± 15 


272 ± 19 


1.05 


-94 ± 24 


301 ± 25 



Mean velocity and velocity disper- 
sion profile. Apertures are 1.00 x 
0.21^^ rectangles. 

and stellar mass. The parameter space is explored using 
a Monte Carlo Markov Chain (MCMC) routine, through 
which the posterior PDF is characterized. The stel- 
lar templates used are based on either a Salpeter or a 
Chabrier IMF. For the description of the photometry of 
the lens we use the double tPIEMD model described in 
Sect. 12. 2[ that is consistent with the analyses presented 
in the following Sections. The stellar masses of the two 
components are fitted independently. Results are listed 
in Table [HI t ogether wit h the values previo usly found by 
iGrillo et all ([2009^ ) and lAuger et al.l (j2QQ9[ ). 

The analysis reveals the presence of dust for compo- 
nent 1, coherently with our previous findings. Repeating 
the fit with the dust-corrected magnitudes yields indis- 
tinguishable stellar masses. The logarithm of the stellar 
masses changes by 0.06 if we use the description of the 
light profile with Sersic components instead of tPIEMDs. 
This is due to the different behavior at large radii of the 
two profiles. Differences in the mass within the outer Ein- 
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Fig. 6. — Mean velocity and velocity dispersion profiles of the 
main lens within 1.15^^ from the centroid. 



TABLE 8 

Stellar mass of the foreground lens, from SPS models 



IMF 



Chabrier 
log(M*/M0) 



Salpeter 
log(M*/M0) 



Comp. 1 
Comp. 2 
Total 



10.85^ 



- + 0.09 
" -0.06 
-I -I 97+0.05 

11.40 ±0.06 



11.13^ 



:, + 0.05 
-0.11 

ii.OZ_Q Qg 

11.66 ±0.06 
11 61+°-°^ 

±±.0±_Q Qg 



This work 
This work 
This work 
Grillo et al. (2009) 



11.34 ±0.12 11.59 ±0.12 Augere^ah (2009 ) 

stein radius for the two models are instead well within 
the measurement errors. 

6. A SINGLE COMPONENT MODEL: MEASURING THE 
AVERAGE SLOPE 
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5200 5300 5400 5500 5600 5700 

Rest frame wavelength (A) 

Fig. 7. — Fit of the velocity dispersion of the lens. Top: The 
red curve is the best fit synthetic spectrum. Shaded regions are 
masked and not used for the fit. Bottom: Residuals of the fit in 
fractions of the total flux. 

In this Section we present a single-component lensing 
and dynamics study of the foreground galaxy, where the 
total density distribution of the lens is described with a 
power-law. The goal is to obtain a measurement of the 
slope of the total mass profile and also to test the accu- 
racy allowed by our data in constraining mass models. 
The system, with its two Einstein rings, offers more con- 
straints than typical single-source lenses. However, the 
analysis is complicated by the presence of two different 
lenses along the line of sight. Light rays from the second 
source are first deflected by the object corresponding to 
the inner ring and then by the foreground lens, with the 
result that, unlike the single lens case, the relation be- 
tween the size of the outer Einstein ring and the enclosed 
projected mass of the lens is nontrivial. Nevertheless, 
this can be properly accounted for as described below. 

A first lens modeling of the system was carried out in 
Paper VI. The procedure adopted there was a conjugate 
points method: multiply imaged spots in the lensed fea- 
tures are identified, and the lens model is determined 
by minimizing the distance between the corresponding 
points in the source plane. This is a conservative ap- 
proach, since it does not make use of all of the informa- 
tion from the surface brightness of the rings. The main 
lens was modeled as a power law ellipsoid, with dimen- 
sionless surface mass density tz = T^/T^cr given by: 



(2) 



where boo = 47r(<JsiE/c)^ and I^is (l^os) is the angular 
diameter distance of the source relative to the lens (ob- 
server). The second lens (first source corresponding to 
the brighter arc) was modeled as a singular isothermal 
sphere (SIS). The model parameter space was explored 
via a MCMC. The results showed that two types of so- 
lution are possible: a model with larger (Jsie, shallower 
slope Y and less massive second lens, or a model with 
a more massive second lens and steeper main lens slope 
(see Figure 9 of Paper VI, or black contours of Figure [H|). 
Part of this degeneracy was due to our ignorance of the 
redshift of the outer ring. 
In this Paper we use the lens model of Paper VI de- 



scribed above and improve it by incorporating 1) our 
measurement of photo-z of the outer ring and 2) a stellar 
dynamics analysis. 

6.1. Stellar dynamics modeling 

We wish to use our measurements of the velocity dis- 
persion profile of the lens to constrain our lens models. 
T his is done with a procedure similar to that adopted 
bv lSuvu et al.l (|2Q1Q[ ), which can be described as follows. 
For a given model provided by the lensing analysis, we 
compute a model velocity dispersion profile and compare 
it to the observed one. The model velocity dispersion is 
obtained by solving the spherical Jeans equation 



1 dp^a^ cr, 



p* dr 



2-^ = 
r 



GM{r) 



(3) 



where p*(r) is the density distribution of the light, and 
(Je are the radial and tangential components of the ve- 
locity dispersion tensor, M(r) is the total mass enclosed 
within the spherical shell of radius r. We impose spheri- 
cal symmetry in the mass model by adopting a spheroidal 
mass distribution 

p(r) (X (4) 

with normalization chosen such that the total projected 
mass enclosed within the Einstein radius equals that of 
the corresponding circularized lens model. The light dis- 
tribution is described as the sum of two tPIEMD profiles, 
with the same parametrization described in Section 12.21 
(best-fit parameters are in table [3]). The 3d stellar dis- 
tribution corresponding to the surface brightness profile 
([Tj) used to fit the photometry is 



p{r) = Perl 



1 



1 



ri + 



rt + 



(5) 



with r = X / + ^ z . Here we set the axes ratios 
q^ to one, as we are assuming spherical symmetry. 

We then assume a O s ipkov- Merr it t model for the ve- 
locity dispersion tensor ()Osipkovlll979l : [Merrittlll985f ): 



= 1 



(6) 



where is the anisotropy radius (orbits are radially 
anisotropic beyond Tq). Finally, we simulate the line- 
of-sight velocity dispersion measured in our apertures. 
Rotation is neglected. Although the lens is seen to be 
rotating, its mean velocity is small compared to the ve- 
locity dispersion and should not contribute much to the 
dynamics of the object. The effect of this approximation 
will be discussed further below. 

6.2. Combining the constraints 

The models of the lens are defined by the set of pa- 
rameters rj = {crsiE,iens, crsis,si , ^s2}: the strength and 
power-law index of the foreground lens, the strength of 
the background lens and the redshift of the outer ring, re- 
spectively. Each model gives a prediction of the velocity 

dispersion in each aperture, o-^^^\ The new posterior 
probability distribution for the model is obtained via im- 
portance sampling: the MCMC sample corresponding to 
the lens modeling of Paper VI is weighted by the likeli- 
hood of the measurements d = {^s2, cri^^^^^^} given the 
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model parameters rj. The following likelihood function 
is used: 

L{d\rj)=P^{z,2)l[G,{cT^:;:r^\v) (7) 

i 

where Pz{zs2) is the PDF in Figure H] and 

/ (meas) _ (mod)N 

(8) 

and crap,i and Acr,i are the zeroth and second moment of 
the posterior PDF of the measured velocity dispersion in 
aperture z, respectively. 

In Fig. [8]we show the updated Posterior PDF obtained 
by importance sampling with the photo-z and dynamics 
measurements, both separately and jointly. It is clear 
that although photo-z and stellar kinematics alone leave 
some degeneracies, the posterior pdfs are almost perpen- 
dicular in this space, and therefore the combination of 
the two is particularly effective. The estimate of the 
slope obtained by marginalizing over the other parame- 
ters is 

y = 1.98 ±0.02. (9) 

We stress that our uncertainty on this parameter is a 
factor of four smaller than the typical error on 7' from 
studies of sin gle-source gravitati onal lenses with SDSS 
spectroscopy (jAuger et al.ll2010"aL see Figu re [T5l) . Com- 
parable precision was reached by Barnabe et al.l ([2011) 
for a sample of lens systems with two dimensional kine- 
matics constraints from integral field spectroscopy. 

In order to better understand the significance of these 
results, we try to quantify the error introduced by our 
simplified model for the stellar dynamics. Two of our 
assumptions are potential sources of bias: spherical sym- 
metry and the non-rotating approximation. The uncer- 
tainty in the mass determination from kinematics data 
is of order 6a'^/a'^ ~ 10%. Biases on the order of this 
uncertainty or smaller are unlikely to bring significant 
changes to the results of our analysis. By consider- 
ing only the velocity dispersion and neglecting rotation, 
we underestimate the mass of the galaxy by a factor 
^ {vrot/o')'^^ which is within 10% in all apertures but 
one. To gauge the importance of this effect we perform 
the following test. We fit the model velocity dispersion 
profiles to the following "effective velocity dispersion": 
(Jeff = \/cr^ + ^rot)' We then apply the same importance 
sampling procedure described above to get a new con- 
straint on the density slope 7^ The new estimate with 
la uncertainty is: 

y = 1.97 ±0.02, (10) 

which is consistent with the original estimate given by 
([9|). On the basis of this result, we can conclude that our 
approximation of non-rotating halo introduces a system- 
atic error of order 0.01 on the inferred value of the slope 

Quantifying the systematics introduced by the spher- 
ical symmetry assumption is n iore c omplicated. In a 
previous work, iBarnabe et al.l (|2011[ ) performed a ro- 
bust dynamical modelling of 12 SLAGS lenses previ- 
ous ly analysed w i th a sp herical Jeans equation approach 
by I Auger et al.l (|2010a[ ). The slopes 7' inferred by 



Lensing + photo-z 




Lensing + dynamics 




1.95 2.00 2.05 2.10 2.15 2.20 2.25 



Lensing + dynamics + photo-z 




Fig. 8. — Posterior PDF of and asi of the updated (filled 
contours) lens model, together with the old model of Paper VI 
(empty contours). The updated model includes only the photo-z 
measurement of the outer ring in the top panel, only the velocity 
dispersion profile of the lens in the middle panel, and both the 
photo-z and velocity dispersion profile in the bottom panel. The 
levels correspond to 68%, 95% and 99.7% enclosed probability. 
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I Auger et al.l (|2QlQa[ ) are consistent with th e more accu- 
rate measurements of iBarnabe et al.l ()2Qlli ). with a bias 
on 7' of 0.05 ±0.04. However, the uncertainty on 7' that 
we achieve in our work is smaher than that and an esti- 
mate of the bias requires additional work. Two distinct 
effects come into play. First, the lens has a non-circular 
projected shape in both its mass and light distribution. 
Second, the galaxy may even have asymmetries along 
the line of sight. The importance of these effects on our 
analysis is quantified in Appendix B. By relaxing the 
assumption of spherical symmetry the additional uncer- 
tainty on the velocity dispersion is about ba^ j ^ 10%. 
It follows that none of our results change appreciably. 

An independent anal ysis of the system was carried out 
by IVegetti et al.l ()2010l ). The method adopted by them 
is more complex than the one used in Paper VI: they 
made use of information from all the pixels of the lensed 
features to reconstruct the source surface brightness as 
a whole. Using data from the inner ring only, they ob- 
tained the following estimate for the density slope: 

y = 2.20 ± 0.03^'*^*^ (11) 

This is a local estimate of the slope 7', obtained by mea- 
suring the magnification of the arc in the radial direction. 
Our measurement is instead an average slope, obtained 
by fitting a single power-law halo to data spanning the 
lens from the center (dynamics) to the outer lensed ring. 
This difference may suggest that the actual mass dis- 
tribution of the lens is diflFerent from a simple power-law 
halo. It is also for this reason that we proceeded to model 
the system with a more complex model. 

7. A TWO-COMPONENT ANALYSIS: DISSECTING 
LUMINOUS AND DARK MATTER 

We perform a two-component lensing and dynamics 
study where the mass distribution is composed of a dark 
matter halo and a bulge of stars. 

7.1. Lensing and dynamics modeling 

We use a power-law ellipsoid for the dark matter, while 
the stars are described with the double tPIEMD model 
found from the photometry analysis. The second lens is 
again modeled as a SIS. The parameters of the stellar 
distribution are fixed to the best-fit values reported in 
Table [3l The global mass-to-light ratio is left as a free 
parameter, but the relative contribution of the two com- 
ponents is fixed according to the results of the stellar 
population synthesis analysis presented in Section |5l For 
a unit F814W-band magnitude, component 1 is measured 
to be a factor of 1.73 (1.77) heavier than component 2 as- 
suming a Salpeter (Chabrier) IMF. In our lensing model, 
the mass-to-light ratio of component 1 is set to be 1.75 
times larger than for component 2. 

We also allow for constant external shear 7ext with 
position angle PAext and constant external convergence 
Kext in the lens plane. Issues related to the external con- 
vergence are discussed below in a dedicated subsection. 
Compared to the lensing study presented in the previous 
section, this model has two additional free parameters: 
the stellar mass M^^ and the external convergence A^ext- 
Given the very tight constraint on the average slope 7' 
from the single component analysis, we expect to be able 
to determine both the slope of the dark matter halo 7dm 



and the stellar mass M^^ with sufficient accuracy. The 
range of values of the slope of the dark matter halo ex- 
plored in this analysis is 1.0 < 7dm < 3.0. 

The technique adopted to fit the model to the lensing 
data is the same used for Paper VI: a conjugate points 
method implemented with a MCMC. The dynamics anal- 
ysis is carried out with a procedure very similar to the one 
described in § 16.11 we solve the spherical Jeans equation 
for our model and obtain a synthetic velocity dispersion 
profile to be compared to the measured one. The (spher- 
ically symmetric) model mass distribution is obtained by 
circularizing the projected mass distribution of the lens 
model, setting (^dm and to one, and by taking the 
corresponding spherical deprojections. The light distri- 
bution is set by circularizing the double tPIEMD profile 
specified in Table [3l 

We then proceed to incorporate information on stellar 
dynamics and on the redshift of the background source. 
This is done by importance sampling, with the same 
method described in § 16.21 

7.2. External convergence 

Objects other than the main lens can contribute to the 
surface mass density k,. This external convergence is hard 
to detect and is degenerate with the total mass of the lens 
galaxy. Ignoring the contribution to from perturbers 
can lead in principle to a bias in the measurement of the 
key parameters of the lens. In order to take into account 
the effect of external convergence on our error budget, 
we include it in our model by generating random values 
of 

/^ext drawn from a plausible distribution. This pro- 
cedure allows us to propagate correctly this uncertainty 
to the other model parameters. Kinematics information 
can also help to constrain A^ext to some extent, as it is 
only sensitive to the mass dynamically associated with 
the galaxy, in contrast to lensing that is sensitive to all 
mass structures along the line of sight to the source. 

Insight on the actual value of /^ext can be gained 
by studying; the lens environment. According to 
iTreu et aTl (12009), this is found to be marginally un- 
derdense with respect to average lines of sight, there- 
fore there is no evidence for the presence of a group 
in the lens neighborhood. The closest cluster known 
to the NASA/IPAC Extragalactic Database (NED) is 
MaxBCGJ146.87912+10.07800, at redshift z = 0.151 
and projected d istance 8.70 arcmin from our lens 
(jTreu et al.l l2QQ9f ) . If we assume a SIS profile for the 
cluster with a typical value for its velocity dispersion 
a = 1000 km s~^ we obtain a contribution to the conver- 
gence /^ci < 0.01. We also scanned the Sloan Digital Sky 
Survey archive looking for massive red galaxies within 5' 
of the lens. Only one early-type galaxy was found, at a 
redshift z = 0.218 and angular distance 2.6'. If we as- 
sume that this object is the brightest galaxy of a group 
and associate it with a SIS halo of <j = 500 km s~^ the 
corresponding convergence at the locatio n of the lens is 
K = .02. Finally, the lensing analysis of iGavazzi et al.l 
(j2008f ) quantified the external shear as 7ext = 0.07 di- 
rected —31 degrees East of North. The HST images show 
two objects with the same alignment relative to the lens 
(see Fig. [3]). If we make the assumption that those ob- 
jects are responsible for the shear and assume again a 
SIS profile we obtain /^ext = |7ext| = 0.07. 
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iHilbert et al.l (|2QQ7[ ) studied the external convergence 
associated with strong lensing systems in cosmological 
simulations. They found that for a source at redshift 
Zs = 5.7 the distribution of A^ext is skewed with a peak 
at —0.04, has zero mean and a scatter of 0.05. A 
slightly smaller sc atter and a peak at —0.02 is found by 
iSuvu et all (j20Tol ) for sources at Zg = 1.39. 

Taking all these aspects into account, we adopt as prior 
for /^ext in OUT analysis a Gaussian distribution peaked at 
0.05, with dispersion = 0.05 and truncated to values 
in the interval —0.05 < A^ext < 0.15. This range should 
capture the indication of a positive contribution from the 
object responsible for the shear and take into account 
the effect of random mass clumps along the line of sight. 
Priors with a broader range of allowed values of /^ext lead 
to larger uncertainties on the other model parameters, 
but none of the conclusions of our study is altered. 

7.3. Results 

Contour plots of the posterior PDF for the model pa- 
rameters are shown in Figures [9] and [TOl The best-fit 
velocity dispersion profile is plotted in Figure [121 The 
inference on the two key parameters and jdm is 
shown in better detail in Figure [Til By marginalazing 
over the remaining parameters, our model constrains the 
stellar mass to 



TABLE 9 

Stellar mass of the foreground galaxy 



M* =5.5t? j X IO^^Mq. 



(12) 



This estimate comes from lensing and dynamics data, 
and does not rely on assumptions on the mass-to-light 
ratio of the stars. This value will be compared with the 
measurement of the stellar mass obtained independently 
from photometry. 

Another important result is the constraint that we ob- 
tain on the slope of the dark matter halo: 



7DM = 1.7 ±0.2 



(13) 



This result shows strong evidence for a contraction of the 
dark matter distribution relative to the inner slopes 
typic al of dark matter only simulations (jNavarro et al.l 
I1997D . Figure [HI shows the mean density profile of each 
mass component compared to the mean single power-law 
fit from Sect. [6l 

Our inference for the anisotropy radius constrains 
> 13 kpc, meaning that radial anisotropy is ruled 
out in the region probed by our data. This is consis- 
tent with pre vious work (e.g. Koopmans fc Treul 120031 : 
iTreu fc Koop mans 2004) and expected on theoretical 
grounds, because strong radial anistropy would lead to 
instabilities. 

In contrast to the power-law model considered in the 
previous section, this mass model has a density distribu- 
tion with a slope that changes with radius. It is interest- 
ing to compare the local value of the slope at the location 
of the i nner rin^ with themeasurement of Vegetti et all 
(|20Ta . iVegetti et all ([20Tol ) modeled the HST F814W 
image using only lensing information from the inner ring. 
Lensing is only sensitive to projected masses, therefore, 
in order for the comparison to be meaningful, we have 
to consider the logarithmic slope of the total projected 
mass distribution, evaluated at the inner Einstein radius. 
We find 

§;^ = -l.l±0.1. (14) 

alogr 



Method 


M* (Mq) 




Lensing+dynamics 






SPS, Chabrier IMF 


(2.5 ±0.3) X lO^i 


2.0 ±0.4 


SPS, Salpeter IMF 


(4.5 ±0.6) X lO^i 


1.1 ±0.2 



^ q; is the IMF mismatch parameter defined as o; = 

This value i s con sistent with the slope found by 
IVegetti et all (j2010D . which is given by -(7' - 1) = -1.2, 
where y is the slope of the 3d mass distribution given in 

m- 

Finally, it is interesting to note how the inference on 
the stellar mass is rather insensitive to the actual value 
of the redshift of the second source, Zs2 (see Figure [9]). 
This means that, with the current data quality, a spec- 
troscopic measurement of the redshift of the outer ring 
would not bring significantly more information. 

8. DISCUSSION 
8.1. Luminous and dark matter in the lens 

The data in our possession allowed us to study the 
lens galaxy of the system SDSSJ0946+1006 under mul- 
tiple aspects. Thanks to the high resolution photometry 
from HST we were able to note now the light distribution 
is well described with two components, while single com- 
ponent models yield poor fits. These two components 
appear to be nearly perpendicular (in projection), have 
significantly different effective radii and surface bright- 
nesses. The colors of the more compact component (com- 
ponent 1 from now on) are also significantly redder (see 
Table Hj), indicator of an older or more metal-rich stellar 
population. As we will discuss below, these characteris- 
tics suggest a particular scenario for the past evolution 
of this object. 

In Sections [Hand [71 we presented two independent mea- 
surements of the stellar mass of the foreground galaxy 
of the system SDSSJ0946+1006 derived with a lens- 
ing+dynamics analysis and with a stellar population syn- 
thesis study. The measured values of M*, obtained by 
marginalizing over the other model parameters, are re- 
ported in Tabled 

The stellar mass measured from gravitational lens- 
ing and dynamics, M^^, is larger than the masses ob- 
tained from the SPS study, Mf^^. This discrepancy 
can be quantified with the "IMF mismatch" parameter 
a = Ml^ /M^^^, also reported in Table H A Salpeter 
IMF is clearly favored, while the probability of the IMF 
being heavier than Chabrier (achab > 1) is 95%. 

This res ult is in agre e ment w ith a gener a l tren d ob- 
served bv iGrillo et al.l (j2QQ9[ ). iTreu et al.l (j2010D and 
lAuger et al.l (|2010bl ) for the early-type galaxies of the 
SLAGS sample. They find that, on average, a Salpeter 
IMF better matches the measurements of stellar masses 
from lensi ng and dynamics. A similar result is found 
by Spinie llo et al.l (2011) for a very massive ea r ly-typ e 
galaxy. As discussed extensively by iTreu et al] (j201Q[ ) , 
stellar mass and slope of the dark matter halo are degen- 
erate with respect to typical lensing and dynamics con- 
straints: given a bulge-halo decomposition, steepening 
the dark matter profile and decreasing the stellar mass 
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Fig. 9. — Posterior PDF in the multidimensional space spanned by the stellar mass M^^, slope of the dark matter halo 7dm 5 radial 
anisotropy scale radius ra, strength of the second lens asi and redshift of the second source Zs2- The levels correspond to 68%, 95%, 99.7% 
enclosed probability. Solid contours: constraints from lensing only. Shaded regions: constraints from lensing, dynamics, and photo-z. 

can result in fits to the observed velocity dispersion and 
mass within the Einstein radius as good as the original 
model. ITreu et al.l (|2010) explained how the observed 
trend of increasing a with velocity dispersion can either 
be interpreted as the effect of a correlation between IMF 
or dark matter inner slope with total mass. Auger et all 
(j2QlQb[ ) explored this de generacy by considering adia- 
batically contracted DM halos set by an imposed rela- 
tion between stellar and virial mass, and found prefer- 
ence for a stellar mass-to-light ratio closer to a S a lpeter 
than a Chabrier IMF. Similarlv. iNapolitano eFall (|2Q1Q[ ) 
find that a Kroupa IMF, which has a mass-to-light ratio 
slightly larger than a Chabrier IMF, fits well adiabati- 
cally contracted DM halos. In the present study we al- 
lowed the slope of the dark matter halo of our lens galaxy 
to vary freely. Its measured value, 7dm = l-7lo!25 sig- 
nificantly steeper than the inner slope of a NFW halo. 
Still, we find a stellar mass larger than what can be ac- 
counted for with a Salpeter IMF and not compatible with 
a Chabrier IMF. Our results imply that a Salpeter IMF 
provides a far better description of the mass-to-light ra- 
tio of the stellar population than a Chabrier IMF even 
with a steepened dark matter h alo. This result is consis- 
tent with the recent findings of iCappellari eFall f2Q12 ) 



and van Dokkum fc Conrovl (|2QTT1 ). 

In contrast, Salpeter-like IMFs are typically r uled out 
for lower mass systems (Cappellari et al.ll2QQ6D or spi- 
ral galaxies (Bell fc de Jong 2001: iDutton et all 120111 : 
ISuvu et al.ll2011l : iBrewer et al.ll2012l ). 

The lensing and dynamics analysis presented in Sect. 
[71 showed evidence for contraction of the dark mat- 
ter halo with respect to a baryonl ess N FW profile. 
A similar result is found by Grillol (120121) for an en- 
semble measurement of 39 massive elliptical galaxy 
lenses. This result is in qualitative agreement with 
many the oretical studies of the evolution of spheroidal 
galaxies (Blumenthal et al." '1986"; 'Gnedin et al.' '20M 
Gustafsson et al. 2006; Abadi et al. 2010; Duffv etldl 
2010). Duffy et al.j (|2010D in their simulations of red- 
shift z = 2 galaxies find inner dark matter slopes that 
span the range 1.4 < 7dm < 2.0 depending on the dif- 
ferent prescriptions adopted to model the effect of the 
baryons. Our measured value of 7dm falls nicely in that 
range, although our galaxy is at significantly lower red- 
shift. Gnedin et al. (2004) provide a prescription to cal- 
culate the dark matter profile of their modified adiabatic 
contraction (MAC) model. It is interesting to test the 
MAC model on the measured slope of the dark matter 
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Fig. 10. — Posterior PDF in the multidimensional space spanned by external convergence ^Cext, strength and position angle of the external 
shear, 7ext, PAext, axis ratio of the dark matter halo qdm, position angle of the major axis of the dark matter halo, PAdm- The levels 
correspond to 68%, 95%, 99.7% enclosed probability. Solid contours: constraints from lensing only. Shaded regions: constraints from 
lensing, dynamics, and photo-z. 
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Fig. 11. — Posterior PDF projected in the space M* — 7dm- 
The vertical shaded regions show independent measurements of 
the stellar mass from photometry, presented in Section 5. 



halo of our galaxy. The final dark matter density profile 
of the MAC model of .Gnedin et ah. (2004) is determined 
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Fig. 12. — Best-fit velocity dispersion profile of the lens. Solid 
line: two components model. Dashed line: single power-law model. 

given the observed light profile, the concentration param- 
eter c of the original (non contracted) NFW halo and the 
baryon mass fraction within its virial radius, f^. 



14 



Sonnenfeld et al. 



10^ 



10^ 



Stars+DM 

DM 

Stars 

Single component 



icr 



10" 
r (kpc) 



1012 



lOii 
>10i" 



10^ 



10« 



10 




Stars+DM 

DM 

Stars 

Single component 



10" 
r (kpc) 



Fig. 13. — Best-fit density (top) and mass (bottom) profiles. 
Solid line: total mass from bulge-halo decomposition. Dashed 
line: stellar mass. Dotted line: dark matter. Dash-dotted line: 
total mass from single component analysis. The shaded regions 
represent 1 — a uncertainties. 

Since we do not have information on the initial prop- 
erties of the dark matter halo of our galaxy, we use a few 
trial values of the virial mass Mvir, spanning a plausi- 
ble range indicated b y a weak lensing study of ellipticals 
(|Gavazzi et alJl2QQ7f), and emplo y a mass-concentration 
relation from lMaccio et all (j2QQ8f ) based on WMAP5 cos- 
mological parameters. We then calculate the inner slope 
of the final dark matte r di stribution with the software 
Contra (Gnedin et al.j l2QQ4f ). The inferred inner slope 
for log (Mvir/Mo) = 12.0, 13.0, 14.0 is plotted in Fig. M 
Despite the large range of virial mass explored, the slopes 
of the contracted halos lie around 1.5 < 7dm < 2.0 over 
the spatial region covered by our data. The MAC model 
is therefore able to reproduce our measurement of the 
dark matter halo slope. 

8.2. A formation scenario 

As our data show, the stellar distribution in the lens 
galaxy consists of two components that differ in align- 
ment, surface brightness and stellar population. This 
particular structure suggests different formation histories 
for the two components. The bright and compact compo- 
nent may have formed first, and later on accreted stellar 



3.0 



2.5 



a 2.0 



1.5 



Inner slope of dark matter halo from MAC 




Fig. 14. — Solid lines: Inner slope of the dark matter halo for 
modified adiabatic contraction (MAC; Gnedin et al. 2004) models. 
Dashed lines: Slope of the non-contracted (NFW) dark matter 
halo. Shaded region: 68% confidence interval of the slope measured 
in this Paper. 

systems in the outskirts without disrupting the structure 
of the original bulge. Alternatively, component 2 might 
have been present originally and component 1 be formed 
in a star formation event following a wet merger. We 
point out that in the infrared image we see evidence for 
tidal distortion in the outskirts of the galaxy (see Fig- 
ure [3]), possible indication of an ongoing merger. Part 
of the faint extended envelope of component 2 could be 
material accreted relatively recently. The presence of 
the dust lane in the center of the galaxy (see § [A]) may 
also be the re sult of a recent merger. We also note that 
IVegetti et al.l (pOlO) detected a compact substructure of 
mass ~ 3 X lO^M© located in the proximity of the inner 
ring image, indicating that minor mergers may still be 
occurring. 

Let us consider our first hypothesis: the galaxy con- 
sisted initially of the compact component 1. What are 
the structural parameters of component 1 and how does 
it relate to other elliptical galaxies? Its effective radius 
is Teff = 0.50'^ (see Tables [2] and [3j) , which corresponds to 
a physical radius of 1.79 kpc. Similar effect ive radii are 
found for hi^h redshift (z > 1.2) ellipticals (iDaddi et al.l 
120051 : iTruiillo et al.l 120061 : Ivan Dok kum et al.ll2008[ ). Its 
stellar mass as inferred from the SPS analysis is given by 
log(MSPS/M0) = 10.85 (log(MSPS/M0) = 11.13) for 
a Chabrier (Salpeter) IMF. Local ellipticals with simi- 
lar values of the stellar mass have effe ctive radii a fac- 
tor of a fe w larger than this object (jShen et al.l 120031 : 
iHvde fc"B ernardi 2009). Analog ously, the high red - 
shift objects of Daddi et al. (20051). ITruiillo et al.l (j2006f ) 
and van Dokkum et al . (2008) are also significantly more 
massive than local galaxies with similar effective radii. 
Finding objects in the local universe that correspond to 
these high redshift "red nuggets" is in fact a standing 
problem in the study of elliptical galaxies. It is not clear 
how objects initially so compact evolve into the more 
diffuse galaxies that we observe at recent times. 

Recent numer ical simulations (|Hopkins et al.l I2009al : 
lOser et ani2012f ) showed how minor dry mergers can in- 
crease the size of elliptical galaxies significantly, with the 
stars of the accreted objects that grow the outskirts of the 
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galaxy, even though the observed and predicted merger 
rates are such that this n iechanism might not be suffi- 
cient ewman et alJl2Q12f ). The observational signature 
of this process would be the presence of a compact core, 
the original red nugget, surrounded by a more diffuse dis- 
tribution of stars from the accreted systems. The galaxy 
studied in this paper might be one of these objects. 

9. SUMMARY 

We have presented a new set of photometric and 
spectroscopic data for the gravitational lens system 
SDSSJ0946+1006. We used these data to constrain the 
structural properties of the foreground elliptical galaxy 
of the system. On the basis of our results, the following 
statements can be made. 

• The redshift of the source corresponding to the 



• The particular structure of the stellar distribution, 
with a comuact core and a misalie:ned faint ex- 



outer ring is Zs2 = 2.41 



+0.04 
-0.21 



at 68% confidence 



level, as revealed by our photo-z measurement. 

• If we describe the total mass distribution with a 
power-law ellipsoid p oc , lensing and dynamics 
data give as measured value 7' = 1.98 ±0.02 ±0.01. 
This parameter should be interpreted as an effec- 
tive slope of the density profile averaged over the 
region within the outer Einstein ring. The spe- 
cial lensing configuration and the exquisite data 
quality of our data, allowed us to measure 7' 
with unprecedented precision. The value obtained 
is consistent with isothermal (7' = 2) and is in 
agreement with the general trend observed for the 
massive early- type galaxies of the SLAGS sample, 
(7') = 2.078 ± 0.027 with intrin sic scatter cr^/ = 
0.16 ± 0.02 (lAuger et all 1201 Oat iKoopmans et al.l 
[2OOI iBarnabe et al.l 1201 ih . See Figure [T5] for a 
comparison of our measurement of 7' with measure- 
ments of the same param eter for the SL AGS sample 
of early- type galaxies bv lAuger et al.l (|2010al) . 

• We are able to decompose dark and stellar matter 
with lensing and dynamics data, assuming a power- 
law density profile for the dark matter. The derived 
stellar mass is S.Sl^ g x lO^^M©, consistent with 
a Salpeter IMF and inconsistent with a Ghabrier 
IMF. This constraint on the IMF is plotted in Fig- 
ure [m together with similar measurements for the 
other SLAGS lenses obtained by Treu et al. (2010). 
Note that we achieve better precision despite using 
less strict assumptions on the dark matter profile. 

• The slope of the dark matter halo is found to be 
7dm = 1.7 ±0.2. This is a strong evidence for con- 
traction relative to the behavior of NEW pro- 
file observed in simulations without baryons, and 
is in agreement with the inn er dark matter pro- 
files obtained by iDuffv et al.l (|201 0) in their simu- 
lat ions of z = 2 galaxie s and with the MAG model 
of iGnedin et al.l (|2004f ). Our inferred bulge-halo 
decomposition has a local projected slope at the 
in ner ring in agreemen t with the value measured 
by "Vegetti et al." ('2010') based on a completely in- 
dependent technique. 
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Fig. 15. — Top panel: IMF mismatch parameter a = M^°/Mf 
relative to a Salpeter IMF vs. lens strength crsiE for the SLACS 
lenses of Treu et al. ( 2010) (black crosses) and for the Jackpot (red 
cross). Measurements of Treu et al. (2010) are obtained assuming 
a NFW dark matter halo with fixed scale radius for the lensing 
and dynamics analysis. Bottom panel: average slope of the total 
density profile 7^ v s. lens strength crgiE for the SLACS lenses of 
lAuger et al.l (|2010al ) (black crosses) and for the Jackpot (red cross). 

tended envelope, might be the result of accretion 
of low mass systems by a compact red nugget. 

• A spectroscopic detection of the redshift of the 
outer ring would still help improve the model, but 
would not lead to a dramatic change in the results 
of our analysis. 
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Fig. 16. — Dust correction. Left panel: Dust map of the central part of the galaxy based on F438W, F606W and F814W images. Dark 
pixels represent higher column densities of dust. Middle panel: Original image in the F336W band. Right panel: Dust-corrected F336W 
band image, showing a single clump of light. The central parts have lower signal-to-noise ratio as a result of the noisy dust map. 



APPENDIX 
DUST CORRECTION 

The presenc e of dust comp li cates our a nalysis. Neverthel ess, we deal with it by applying a procedure similar to that 
adopted by Koop mans et al.l (|2QQ3[ ) and iSuvu et al.l (|2QQ9[ ) for the system B1608+656. The details of the procedure 
are the following. 

We select a small region in the center of the galaxy for which we want to apply a dust correction. We estimate the 
intrinsic colors of this central p art by measur i ng the m in a region that we think is free of dust. 

We assume a dust law from iCardelli et all (|1989f ). with Ry = 3.1. Given the flux in one band and the colors of 
the object, the fluxes in the remaining bands are determined by the dust law. Therefore, with images in two or more 
bands we can constrain both the intrinsic flux and the dust content of the object. In our case we determine these two 
quantities in the central region of the lens on a pixel-by-pixel basis by fitting the F814W, F606W and F438W fluxes. 
We want to check if we can account for the dark lane observed in the F336W image independently from the data in 
that band, therefore we do not include that image in our dust analysis. The F160W image is excluded because of its 
lower spatial resolution. The PSF of the different images is not matched. Our inference on the presence of dust is not 
affected by this approximation. 

The two dust-free colors, F606W - F814W and F438 - F606W, are measured in an annulus around the center and 
inside of the inner ring. We cannot rule out the presence of a uniform dust screen, but that would not affect our 
conclusions as the tools that we use for quantitative analyses can account fo r that. The d ust-co rrected flux in the lens 
center is then calculated with the fitting method described in Sect. 5.4 of iSuvu et al.l (|2QQ9f ). Figure [16] shows the 
recovered dust map, the F336W image corrected for dust and its original version. It can clearly be seen how the dust 
map, obtained without using data from the F336W, has largest column density right where we observe the dark lane 
in the image. In the dust-corrected image, the lens looks indeed more like a single object. The overall amount of dust 
is relatively small, as the correction to the total magnitude in the F606W band is approximately 0.10 mags. 

ELLIPTICITY EFFECTS 

The mass and light distributions of the lens galaxy are well approximated by ellipses, but we make the assumption 
of spherical symmetry for the analysis of the kinematics data. How does this assumption a ffect the results we presen t? 
To answer this question we make use of the axisymmetric version of the Jeans equations (|Binnev fc Tremaind 120081 ) 

v^{R,z)-v^ + - +R—, (B2) 

obtained assuming a distribution of stars of the form f{E^ L ^). 

We take our best-fit single component model from Section 16.11 and make it elliptical by assuming that the rotation 
axis Lz is in the plane of the sky and fixing the projected ellipticities in light and mass to = 0.95 and q = 0.87 
respecti vely. The first value is the ellipticity of the best single-component fit to the light profile, the latter value is 
given in iGavazzi et al.l ()2008f ). Then we assume isotropy in the velocity dispersion tensor, calculate the line of sight 
velocity dispersion profile in the two possible cases of oblate or prolate ellipsoid, and compare it to the corresponding 
spherical case. Results are shown in Fig. [Bl Deviations from spherical symmetry bring differences on the order of a 
few km s~^ on the velocity dispersion profile, well within our uncertainties on the measurements, and therefore is not 
a concern for possible biases. 

More important are the effects of asymmetries along the line of sight. We do not have any direct measurement 
of the line-of-sight structure of the lens, but from t he observed proje c ted fla ttening we can get information on the 
intrinsic shape of the galaxy by statistical means. iPadilla fc StraussI (|2008 h measured the distribution of intrinsic 
axis ratios of massive elliptical galaxies. By drawing samples of galaxy shapes from their inferred distribution and 
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Fig. 17. — Top: Line-of-sight velocity dispersion profile for a spherical model and for elliptical models with = 0.95 and q = 0.87, 
calculated by solving the axisymmetric Jeans equation. The models have the same projected mass within the inner Einstein radius. PSF 
smearing is not included, resulting in the high central peak. Overplotted are the error bars on the measured velocity dispersion profile. 
Bottom: Velocity dispersion profile of the spherical model and of an oblate (prolate) ellipsoid with minor (major) axis parallel to the line 
of sight and axis ratio 0.8. 

assuming random orientations we find that 68% of the objects that produce a projected ellipticity = 0.95 have 
an axis ratio rounder than 0.8. How does the velocity dispersion profile of an oblate (prolate) galaxy with minor 
(major) axis along the line of sight and axis ratio of 0.8 differ from that of a spherical galaxy with the same (observed) 
projected mass within the Einstein radius? We use the axisymmetric Jeans equation to address this question as well. 
We take our best-fit spherical model and modify it into an oblate (prolate) ellipsoid with the axis ratio of both the 
light and mass distribution fixed at 0.8, orienting along the line of sight. The line-of-sight velocity dispersion profile 
for isotropic orbits in the oblate and prolate case is also plotted in Fig. [B1 The spread relative to the spherical case 
is somewhat larger than the uncertainties. To make sure that our assumption of spherical symmetry does not alter 
the measurements presented in this Paper we recalculate the inference of the key model parameters by inflating the 
error bars on the velocity dispersion measurements by a factor 1.5, matching the scatter introduced by the unknown 
line-of-sight oblateness or prolateness of the lens. None of the results changes appreciably. The lens modeling does 
not depend on assumptions on the line-of-sight mass distribution, and so in this regard our results are robust. 



